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Abstract 

In this paper, stability of linearly coupled dynamical systems with feedback pinning algorithm is studied. Here, both 
the coupling matrix and the set of pinned-nodes vary with time, induced by a continuous-time Markov chain with 
finite states. Event-triggered rules are employed on both diffusion coupling and feedback pinning terms, which can 
efficiently reduce the computation load, as well as communication load in some cases and be realized by the latest 
observations of the state information of its local neighborhood and the target trajectory. The next observation is 
triggered by certain criterion (event) based on these state information as well. Two scenarios are considered: the 
continuous monitoring, that each node observes the state information of its neighborhood and target (if pinned) in 
an instantaneous way, to determine the next triggering event time, and the discrete monitoring, that each node needs 
only to observe the state information at the last event time and predict the next triggering-event time. In both cases, 
we present several event-triggering rules and prove that if the conditions that the coupled system with persistent 
coupling and control can be stabilized are satisfied, then these event-trigger strategies can stabilize the system, and 
Zeno behaviors are excluded in some cases. Numerical examples are presented to illustrate the theoretical results. 
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1. Introduction 


Control and synchronization of large-scale dynamical systems have attracted increasing interests over the last 
several decades Jlj]-@|. When a networked system is unstable by itself, many control strategies are designed to 
stabilize the networked system. Among them, pinning control is effective, because it is economically realized by 
controlling a partial of the nodes, instead of all nodes in the network. The general idea behind pinning control is 
that when applying some local feedback controllers only to a fraction of nodes, the rest of nodes can be propagated 
through the network interactions among nodes 0-0. 

More related to the present paper, for example in I18n . the authors investigated the pinning control problem of 
coupled dynamical systems with Markovian switching couplings and Markovian switching controller-set. In 0 and 
most existing works in linearly coupled dynamical systems, each node needs togather information of its own state 
and neighbors states and update them continuously or in a fixed sampling rate Q6D. However, as pointed out in 111911 . 
the event-based sampling technique showed better performance than sampling periodically in time for some simple 
systems. Hence, a number of researchers suggested that the event-based control algorithms can be utilized for the 
purpose to reduce communication and computation load in networked systems [{20;]- [ 21 ] but still maintain control 
performance [211- 12711 . Therefore, the event-based control is particularly suitable for networked systems with limited 
resources and so has attracted wide interests in the scope of distributed control of networked systems. The idea 


‘Corresponding author. 

Email addresses: wenlianS fudan. edu. cn (Wenlian Lu), y jhanSshmtu .edu.cn (Yujuan Han), tchen@ fudan. edu. cn 

(Tianping Chen) 


Preprint submitted to Elsevier 


February 3, 2015 










of event-triggered control can be regarded as a specific temporal discretization approach, which was studied before 
Hi. As application, the event-based control was applied for consensus of multi-agent systems. For instance, 
investigated centralized and distributed formulation of event-driven strategies for consensus of multi-agent systems 
and proposed a self-triggered setup; [23, 3 studied the stochastic event-driven strategies; [25] introduced event- 


based control strategies for both networks of single-integrators with time-delay in communication links and networks 
of double-integrators; By using scattering transformation, l26ll studied the output synchronization problem of multi¬ 
agent systems with event-driven communication in the presence of constant communication delays. 

In some recent papers [2 ( )|-|3()|. the authors addressed event-triggered algorithms for pinning control of networks. 

| gave an exponentially decreasing threshold function, while the event-triggering threshold in 128. 3(ih is prescribed 
by the continuous or discrete states of agents and target. In these works, sufficient conditions were proposed, which 
are based on the control gain, some quantities of the uncoupled node dynamics and the minimum eigenvalue of the 
augmented Laplacian. 

Motivated by these works including our previous work & in the present paper, we employ the event-triggered 
strategy in both coupling configuration and pinning control terms to realize stability in coupled dynamical systems 
with Markovian switching couplings and pinned node set. At each node, the diffusion coupling and feedback pinning 
terms are piecewise static based on the information of its local neighborhood and the target trajectory only at the latest 
time of event, which is triggered by some specified criteria derived from the information of its local neighborhood 
and target. In other words, once the triggering criterion of node is satisfied, the diffusion coupling and pinning 
terms will be updated; otherwise these terms are constant between two successive event time points. We consider 
two scenarios: continuous monitoring and discrete monitoring. In the continuous-monitoring scenario, each node 
observes its neighborhood’s and the target’s states in an instantaneous way; on the contrary, in discrete-monitoring 
scenario, each node can only obtain its neighborhood’s and target’s information at the last event-triggering time point, 
which results in a small cost of monitoring (communication load) but the triggering events happening more frequently 
than continuous monitoring, namely, higher computation load. For both scenarios, it is proved that the proposed 
event-triggered rules guarantee the stability of the coupled dynamical systems under the local pinning algorithm. 

This paper is organized as follows. After formulating the underlying problem in Sec. 2, the event-triggering 
rules of diffusion and pinning terms are proposed to pin the coupled systems to a homogenous preassigned trajectory 
of the uncoupled node system by using continuous and discrete monitoring scenarios in Sec. 3 and 4 respectively. 
Simulations are given in Sec. 5 to verify the theoretical results. Strength, limitations of the work and possible orients 
of future study are discussed in Sec. 6. Finally, this paper is concluded in Sec. 7. 


2. Problem formation 


In this paper, we consider a network of linearly coupled dynamical systems with discontinuous diffusions and 
feedback pinning terms as follows: 


Xi(t) 




i I5 * * * i 


( 1 ) 


where 


Oi(vt,t l k ) = -c^L.K)r 

j =1 




- ceDi(a t ) 




( 2 ) 


Here, the system contains to nodes, Xi(t) = [ x\{t ), • • • ,x?(t)\ T £ R" denotes the state vector of node i, the 
continuous map /(•) : E” —> E" denotes the identical node dynamics. a f is a homogeneous Markov chain, which 
will be specified later. L(a t ) = [Lij(crt)]™j=i £ R m,rra is the time-varying Laplacian matrix of the underlying time- 
varying bi-graph = {V, -E(crt)}, with the node set V and time-varying link set F?(<r t ): for each pair of nodes 

i 7 ^ j, Lij(at) = —1 if z is linked to j at time t , otherwise (cr*) = 0, and La(at) = — Lij(&t )■ c is uniform 

coupling strength at each node. F = [ 7 tj]”j =1 £ stands for the inner configuration matrix that describes the 
coupling of components between the state vectors. Di{crt ) = &v(< 7 t ){i)' where <$.(•) is the characteristic function, i.e. 

= 1 if t £ T>(<Jt), otherwise <5x>(o- t )(*) = 0 for the pinned node subset T>(a t ) C {1, • • • , to}, where T>(a t ) 
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denotes the node subset in V that are pinned at time t by a specific node dynamic trajectory s(t) with s = f(s(t)), 
s(0) = So- e is the pinning strength gain over the coupling strength. 

Our aim is to provide sufficient conditions to guarantee that s(t) is a global stable trajectory for the coupled 
system, namely 


lim ||(t) — s(t) || = 0, V i = 1, • • • , m. 

t—y OO 

Here, we consider the /^-vector norm and denote it by || • || throughout this paper. Let Xi(t) = Xi(t ) — s(f). Then C] 
and |2]) become 


Xi(t) = f(xi(t) + s(t)) - f(s{t)) - cy^ j L ij (a t )r 


j =i 




- ceDi(a t )Txi(tl), t\<t< t\ +1 (3) 


Suppose at is a homogeneous continuous Markov chain with a finite state space S = {1,2, ■■ ■ , N } and its 
infinitesimal generator Q = [q uv ]NxN is given by 


P{ct*+a = v\a t = u} 


q uv A + o(A), u^v, 
1 + q uu A + 0 (A), u = v, 


where A > 0, o(A) is a infinitesimal as A — > 0, i.e, limA->.c)(o( A )/ A ) = 0, p uv = — > 0 is the transition 

probability from u to v if v u, while q nu = — Y^v=i v =£u Denote P = [p uv \ the transition matrix of the 
Markov chain. The sojourn time in state u is exponentially distributed with parameter q u = —q uu . 

Let Pt ill (s, A) = P(tJt+ s £ A|oy = u) and E t . u (/(cr s )) = f f(y)Pt,u(s,dy). Denote by A the weak infinitesimal 
operator of a t . A function /(<r t , t) is said to be in the domain of A if 

lim Et,«(/(g t+ A,£ + A))-/(u,£) 

A-vO A 

E t , u (f(a t+ A,t + A)) - E t ,„(/(cr t+ A,£)) ,. E t , u (/(cr t+A ,£)) -/(u,£) 

A—>o A A-s-0 A 

=ft(u,t ) + h(u,t) 


and 


lim E t „(/ t (CT t+ A,f + A) + h(a t +A,t + A)) = f t (u,t) + h(u,t). 

A—>0 

Then, we write Af(u, t) = ft(u,t ) + h(u, t). 

Note that for fixed t, 

Et,u(/(CT t+ A, £)) = f f(y, t) P t ,„(A, dy) = ^2 f ( v ) v ) = ^ quv A + u > *)■ 

V V 

Hence, by the Dynkin’s formula, we have 

Af{u,t) = ft(u,t) +^2f(v,t)quv (4) 

V 

Throughout the paper, we assume /(•,£) belongs to the following function class. 

Definition 1. Function class QUAD (G, ak, 3): let G be an n x n positive definite matrix and T be an n x n matrix. 
QUAD(G, aT, {$) denotes a class of continuous functions /(£, t) : R” x [0, +oo) K > R n satisfying 

(£ - C) T G[/(e, t) - /(C t) - aT(£ - C)] < -Pit ~ C) T G(C - C) 


holds for all £, £ £ R™. 
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Definition 2. System ( 17121 > is said to be exponentially stable at s(t) in mean square sense, if there exists constants 
6 > 0 and M > 0, such that 


E 


\\xi(t) — s(f)|| 2 


< Me~ 6t 


holds for all t > 0 and any * = !,••• , to. 


(5) 


3. Continuous monitoring 


We briefly provide the basic idea of the setup of the coupling and pinning terms. Instead of using the simultaneous 
state from the neighborhood and the target trajectory to realize stability, an economic alternative for the node i is to 
use the neighbors’ constant states at the nearest time point t,' k until some pre-defined event is triggered at time t l k+1 ; 
if node i is pinned at time /, it also obtains the target trajectory’s constant state at time point if (t y, then the incoming 
neighbors’ and the target trajectory’s information is updated by the states at t\ +1 until the next event is triggered, and 
so on. The event is defined based on the neighbors’, the target trajectory’s and its own states with some prescribed 
rule. This process goes on through all nodes in a parallel fashion. To depict the event that triggers the next time point, 
we introduce following Lyapunov function: 

V(x,t,a t ) = ix T ^P(cr t ) ® G^x, (6) 

where x = [x^, • • • , x„] T , P(er t ) £ M m,m are diagonal positive definite matrices, induced by at, and G £ 
is a positive definite matrix. Let F(x) = [(/(xj) - /(s)) T , • • • , (/(x m ) - /(s)) T ] T , D(a t ) = diag[Di(a t )]£L 1; 
L(o t ) = L(a t ) + eD(a t ). 

Note that V (x, t, af) is in the domain of the weak infinitesimal operator of at- Denoting at = u, by ©, we have 


N 


AV(x,t,u ) = ^2 q U vV(x,t,v) + 


dx dt 

V=L 

Substitute ([3]) into i(7}. we get 

AV(x, t, a t ) |3j = x T (t) [P{a t ) ® G] F(x(t)) - al m ®Tx(t) + cz(t) 

+ x T (f) |p(cr t ) (aim - cL{a t )) ® GT + ^ ^ g CTt t,P(x) ® g| x(f), 

where z{t) = [zj ( t ), • • • , zf(t)] T with 

z i(t) = x j (t)-x i (t)-x j {tl i ( t) ) + x i (tl iW ) 


+ eDi(a t ) Xi(t) - s(t) - Xi(t l ki(t) ) - s(t l ki(t) ) 




(7) 

( 8 ) 


(9) 


and denotes the symmetry part of a square matrix S, i.e., S sym = (S + S T )/ 2. 

Let A m (-) and Xm(/) denote the smallest and largest eigenvalues in module of a symmetry real matrix, and 
A = min„ A m (P(v) G), A = max„ A m (P(v) G). From the condition / £ QUAD(P, aT, /?), we have 


x (f) [P(cr t ) g) G] P(x(f)) - aim (8) Tx(t) < — /3Ax (f)x(f). 


For the term x T (f) [ P(at ) ® G\ z(t), we have 

x T (f) [P(a t ) ® G] z(t) < 
< 


^X T (f)(P 2 (cr t ) ® G 2 )x(f) + ^-Z T (t)z(t) 
^-x T (f)x(f) + ^-z T (t)z(t) 


holds for any v > 0. Then, we have the following theorem. 


GO) 

( 11 ) 
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Theorem 1. Suppose that f belongs to QU AD(G, oT, /3) with the positive matrix G and a > 0, /? > 0, and there 
exist diagonal positive definite matrices P(u), u = 1, ■ • • ,N such that 

1 N 

{P(u)[al m — cL(u) — ceD(u)\ (g) GT} sym + - ^ q uv P(v) <S> G < 0, for all u € §. (12) 

V=1 

Then, under either of the following two updating rules, system 0 is exponentially stable at the homogeneous trajec¬ 
tory s(t) in mean square sense: 

(1) set t l k+1 by the rule 

4+1 = max jr > t\ : ||^(r)|| < ■^=^=|^||£’ i (r)||| (13) 

where 0 < S < 2/3X/X is a constant; 

(2) set 4+1 by the rule 


4+1 = max{r > t\ : ||^(r)|| < aexp(-6r)} 
where a > 0 and b > 0 are constants. 


(14) 


Proof [Case (1).] Consider the event-triggering rule lfl3l > and pick a constant S with 0 < 6 < 2/3X/X. By Dynkin 
Formula ifiTtl . we have 

Ee^ V(x, t, a t ) = EV(x(0), 0, oo) + <5E f e Sr V(x, r, o T )dr + E f e ST AV(x,T,a T )dT. (15) 

Jo Jo 

From ©-G3. we have 

E e St V(x,t,<Jt) < EC(a;(0), 0, fJo) + ® f e ST V(x,T,o T )dr — /3AE f e Sr x T (r)i(r)dr 

Jo Jo 

™.\2 ft ft 

H—-—E / e' 5r a: T (r)£’(r)dT 4- — E / e Sr z T (r)z(r)dT 
2 Jo 2 V Jo 

n ~ cL(cr T )) gt + — 

< EV{x{0),0,ao)+E [ e s 

Jo 


f* e Sr x T (t) < 

P(a T ) 

' o 

[ 


sym 


5 X cvX 2 
-+A + —+ — 


x (t)x(t) + — z (t)z(t) > dr (16) 
2v 


for any v > 0. Note that 


2v 

max — 
v>0 c 


SX cvX 2 } _ {pX-\8X) 2 
2 2 


cX 2 


and the maximum is reached if and only if v = ■ Hence, letting v = . (fl3] > implies 


2v 


\zi(r)\\ 2 <- 

c 


y/cX 

SX cvX 2 
~ ~2 2 ~ 


l®i( r )l| 2 ) i = 1, - ■ ■ , TO, 


for all r < 4+i- Therefore, we have 


Ee St V(x, t, a t ) < EC(x(0), 0, <jq), 


(17) 
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which implies 


Ee 5t \\xj(t) - s(i)|| 2 < ^-¥.e St V(x,t,a t ) < yEV(:e(0), 0, a 0 ), Vj = 1, ■ • • ,m. 
A A 


Therefore, we have 


2e 


—St 


E||xj(f) - s(f)|r < —EU(x(0),0,cr o ), Vj = !,-■■ ,m. 


[Case (2).] Consider the event-triggering rule (fl4l) and pick v = ■ Then, we have 

<5A cv A 2 

-P±+ T + — = o. 


Substituting ||2»(t)|| < aexp(— br) into (fl6l) gives 


Me St V(x, t, a t ) < EV(i(0),0,<To)+ 


a 2 c 2 A 2 


= EV(x(0),0,ct o ) ■ 


2(2/7 A — SX) J 0 
a 2 c 2 X 2 1 


2(2/3 A — <5A) 26-5 L 


E / e (- 26+i > T dT 

J _ g( — 2b+S)t 


Let Co = 


0 2(2/3A-<5A) 2b—6 

2e~ st 


i-r. By the similar arguments as m-m, we have 


®j(t) - s(t )|| 2 < (EU(cc(0),0,ct o ) +C 0 [l - e (- 26 +*)*]) < c 1 e- min{2b ’ s '>\ Vj = 1 


(18) 


(19) 


L, , f/lz, 


where C\ = j max(|EV(x(0), 0, op) + Co|, |Co|). This completes the proof. 


□ 


Remark 1. If x(t\) = 0 holds, then from ( I7JD . at time ti, every node updated its feedback term and pinning term (if 
pinned). Therefore, from ( 17121 ). for any t > t\, x(t) = 0 holds, which means x = 0 is an equilibrium of the system 
under the event-triggering rules. 

Remark 2. In fact, in our previous work Mi, we studied pinning dynamic systems of networks with Markovian 
switching couplings and controller-node set. By this theorem, we showed and proved that if the condition for the 
stability of the coupled system with spontaneous coupling and control in M can be satisfied, then the event-trigger 
strategies can stabilize the system too. Therefore, the issue of selection of pinned node in terms of guaranteeing 
stability is totally the same with the system with spontaneous diffusion and control. 

Under the updating rule (IT4l i. it can be proved that the Zeno behaviors 02] is excluded by the arguments as the 
same fashion as in [29]. While for the updating rule (fT3l> . similar to work II2211 . it should be pointed out that there 
exists at least one node such that its next inter-event interval is strictly positive. 

Proposition 1. Suppose that all hypotheses of Theorem\T\hold. Under the updating rule ( I7JI ). if the system does not 
reach stability, then there exists at least one node i € {1, • ■ ■ , to} such that the next inter-event inten’al is strictly 
positive; Under the updating rule ( 1741 ). if the system does not reach stability, the expectation of next inter-event inten’al 
of every node is strictly positive, further, it is lower bounded by some positive constant. 

The proof is the similar to those in [!2§, 22]] with some modifications. In fact, if at any time t, there exists one 
node i such that ( IT3l) cannot hold as an equality for this node i. Hence, t l k+1 > t can be derived, which implies that 
the inter-event interval for node i is positive. Otherwise, at time t. (fT3l> holds as an equality for every node, that is, all 
nodes update their control law at this moment, which implies Zi(t) = 0 holds for all i. However, since the network 
has not been stabilized at s(t) yet, there exists at least node j with Xj(t) ^ 0, which implies t° k+l > t holds, which 
implies that the next inter-event interval for node j should be positive. 

In comparison, under the rule (IT4l >. suppose 6 < |, for each node i, at r = t\, we have Zi(t l k ) = 0 but the right- 
hand side of (IT4l) is nonzero. Therefore, t l k+l always holds and the low bound of the expectation of inter-event 

intervals can be estimated as £ log {1 + A ^ B } with Lf the Lipschtiz constant of /(•), A = 2mL /+ 2cm (™+ <; )+- L /+ cm 
and B = 
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4. Discrete monitoring 


In the discrete monitoring scenario, each node i can obtain its local neighborhood’s state only at the time points 
t\,k = 1,2, • • ■. Meanwhile, if node i is pinned, it can also obtain the target’s state at latest time points t\.y t y By this 
way, the rule to determine the next time point t l k+1 of obtaining state information only depends on the local states at 
i' k . In comparison, the triggering event rules (fl3l) and ( 1 1 41 demand the instantaneous states after t\. 

Consider system (1112b and V(x) as the candidate Lyapunov function with its derivative <[8j. We are to derive a 
triggering event rule from (fT3l> in Theorem [I] which only depends on t\. The estimations of the upper bounds of 
||zi(i)|| and the lower bounds of || Xi(t) — s(f)|| for all i are essential. 

First of all, we take the switching time points of the Markov chain a f to trigger the state information updating for 
all nodes. Then, we are to estimate the upper-bound of \\{xj(t) - Xj(ti)) - ( Xi(t ) - Xi(t \))II with Lijtat) ^ 0, of 
which the evolution equation can be written as: 

/ = /(*<(*)) + W<ru%) r9m 

for t\ < t < min{ij, +1 , t 3 k and initials Xi(t' k ), Xj(t\). Here, 9i(a t , t l k ) and 9j (a t , t k ^) are constants in this 

time interval. 

Let us consider a general form of ( 120b as follows: 


( f = /(«(*))+ 0 «(0) = «o (21) 

1 ^ = + ^ v{0)=v 0 , 

where u, v, uo, vq € K.". Suppose that there exists a nonnegative-valued continuous map p : R>o x R 4 " M>o such 
that the solutions of (fTil satisfy the following inequality: 


H(M(f) - Mo) - (v(t) - Mo)|| < p{t 1 B 1 'd 1 u 0 ,v o). (22) 

Here the map p depends on the node dynamics map /(•), the initial value uq,vq and inputs 0. 0, and satisfies 
p( 0, •, •, •, •) = 0. Geometrically, p is an upper-bound estimation of the difference between the two trajectories of 
(l2ll) starting at 0 of time-length t: 


||(u(f)-u 0 )-(i>(f)-i;o)|| 


[/(u(s)) - f(v(s))]ds + (d- d)t 


For example, if /(•) is Lipschitz (on the two trajectories): ||/(u(s)) — /(m(s))|| < Lf\\u(s) — v(s)|| for all s > 0, 
then we have 


||(u(i) - Mo) - (v(t) - Vq ) 11 <Lf [ II (m(s) - M 0 ) - (v(s) - M 0 )||ds 

Jo 

+ (||6»-r?|| + L f \\u 0 - v 0 \\) t. 

By the Gronwall-Bellman inequality lf33. 341. we have 

ll(«(t) - n 0 ) - (v(t ) - m 0 )|| < (W e -n+L f \\ u o- v o\ l) [exp(£/f) _ !]. 


(23) 


We can take p(t, 9 , •&, Mo, Mo) as the right-hand side above, which equals to zero at t = 0. 

Second, we suppose that there exists a nonnegative map g : R>o x K 4 " —> K>o such that the solutions of (|2TT > 
satisfy: 


\\u{t) - v(t )|| > g(t, 9, M 0 , m 0 ). 


(24) 
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] 



Figure 1: p generator. 


Here, g can be regarded as the lower-bound estimation of the distance between two trajectories: 


IK*) -v(t)\\ 



[f{u (s)) - f(v(s))]ds + (9 - d)t + (u 0 - wo) 


and satisfies (i). g(-, 8 1 6, uq, uq) = 0; (ii). p(0, •, •, ito, ito) = 0. For example, assuming that there exists some 
constant er (possibly negative) such that 


(it - v) T (f(u) - f{v)) > a(u - v) T (u - V) 


holds for all it, v € R". We have 
d 


dt. 


[(u(f) - v(t)) (it(f) - w(f))] J2T} = 2(u - v) [/(it) - f(v) + 8-d] 


1 


> 2cr(it — v) (it — v) — n{u — v) (it — v) - {6 — d) {8 — ■&) 




hold for any g > 0. By Gronwall-Bellman inequality, we have 

(u(t) - v(t)) T (u(t) - v(t)) > exp [(2cr - n)t](u 0 - w 0 ) T (u 0 - v 0 ) 

{8 - d) T (8 - 0)/n 


2a — g 


exp[(2cr — g)t] — 1 >. 


We take g as the right-hand side above, which is positive for a small interval of t, starting from 0, for any uq ^ vq. 

We highlight that there is no uniform approach to get precise estimation for a general function of /(■) but one can 
do it case by case. Therefore, an efficient way is to use integrators that simulates the node dynamics of u = f{u) + 8 to 
realize generators that calculate the maps of p and g. Noting that these generators are independent of the states of the 
nodes, they can be built parallel to the networked system. FiguresQ]and[2]show the generators of p and g respectively. 

Let d l k = 8i(at,t\) and ^ = 8j(a t . Based on the event-triggering rules (fl3l) . (IT4l >. we have the 

following theorem. 


Theorem 2. Suppose that f belongs to QUAD(G , aT, /3) with positive matrix G and a, f3 > 0. Suppose there exist 
diagonal positive definite matrices P(u), u = 1, • ■ ■ ,7V such that 


{P(u){al m 


cL{u) - ceD(u )] ® GT} sym 



N 

q uv P(v) ® G < 0, for all u € S. 

V — l 
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Figure 2: g generator. 


define a sequence of under either of the following two updating rules, 

a) 

Ck = max | f ^ (~ L ij ( a £+t% ))P (f» ^(C+ti )•**(**)’ *j(4)) + 

+eDi ( a i+t i k )p (£,&]., Q,Xi(t\), s(t 2 k )) 

< ^(4)^(4))} (25) 

where 0 < 5 < 2/3A/A is a constant; 

( 2 ) 

Ck = max | £ : ^(-£y(<r e+( j))p (?,)>^(4).^'(4)) + 

+eA(o- e+t i )p (4 tffe, 0, Xi(t l k ),s(tD) 

< a exp (-&(£ +4)) | (26) 


where a > 0 ahc/ b > 0 ore constants. 

If the triggering event time points {t k } are picked by the following scheme: 

1. Initialization: tf = 0 for all i = 1, • ■ • , m; 

2. At t = t\., node i obtains f k by the rule ( 1251 ) for (1261)): 

3. At t > t l k , if one of its neighbor, for example, denoted by j, is triggered at t = t k ,, 1 (^ e! k' be the latest event at 
node j before t), then j broadcasts its current updating law, "Of, to node i, and the rule ( 1251 ) for ( 1261 )) is updated 
by replacing the diffusion term from node j, d k ,, by d k , +1 , and tf by t. And, go to Step 2; 

4. If a t switches at t, then we update the rule ( l25l >ior ( 1261 )) by replacing d\ by the current state 6i(at , t), and t\ 
by t. And go to Step 2; 

5. Let t k+1 = t k + f k , an event is triggered at node i by updating the state information in ( 17 121) from t’ k by t k , 1 , 
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then system ( 17121 ) is stabilized at s(t) in mean square sense. 

This theorem can be derived from Theorem[T]immediately. In fact, event (l25l > is an estimation of event (IT3l) . event 
(l26l) is an estimation of event (IT4l) . 

There is substantial difference between the discrete and continuous monitoring strategies. Generally speaking, the 
continuous monitoring require that every node collects its neighborhood states at every instant time, while discrete 
monitoring does not need this step. As shown in Table Q] the continuous monitoring scheme costs higher commu¬ 
nication load than the discrete monitoring. As a pay-off, we will show in the numerical example section that the 
frequencies of triggering events in the continuous monitoring are much lower than that the discrete monitoring re¬ 
quires. That is, the continuous monitoring costs lower computation load than the discrete monitoring. 


Table 1: Continuous vs discrete time monitoring schemes 


Step 

Continuous monitoring 

Discrete-time monitoring 

1 

At time t\, agent i updates feedback control law 

At time t\, agent i updates feedback control law 

0i(vt,t l k ) 

2 

in (113b or (1 141) 

1F^CI>125Fot126I 

3 

then 

monitoring the states of i’s neighborhood Xj{t), j £ 
Mi and target s(t) (if i is pinned at time t), t > t l k 


4 

else 

go to step 1, replace t k by t\ +1 

else 

go to step 1, replace t k by t\ + f l k 


Remark 3. In the discrete monitoring scenario, each node does not need to observe the information of its neighbors 
at every instants, but each node has to broadcast its updating law, d l k , to all its neighborhood once it is triggered. 

Similar to PropositionQ] we have: 

Proposition 2. (1) Suppose that hypotheses in Theorem\2\hoId. Under the rule ( 1251 ) and the scheme described in 

Theorem\2\ if system ( 17121 ) is not stable at t, there exists at least one node i £ {1, • ■ • , to} such that the next 
triggering event time strictly greater than t, namely, inter-event inten’al is strictly positive. 

(2) Under the rule ( 12(51 ) and the scheme described in Theorem\2\ if system ( 17121 ) does not reach stability, the expec¬ 
tation of next inter-event inten’al of every node is strictly positive, further, it is lower bounded by some positive 
constant. 

Remark 4 . The discrete monitoring strategy implies the triggering events happen more frequently than continuous 
monitoring as a reward of a smaller cost of monitoring. 

Remark 5. For the discrete monitoring strategy, the computation complexity for every task depends on the number 
of multiplies in p(-) and p(-). We suppose the number of multiplies of p(-) and p(-) are respectively N\, N 2 . From 
the updating rules ( 1251 ) and ( 1261 ). the computation complexity for the next triggering time of every agent is at most 
(to + l)7Vi + N 2 . 


5. Examples 


In this section, we present several numerical examples to illustrate these theoretical results. The system is an array 
of 5 coupled Chua circuits with the map /(•) of node dynamics as follows: 


m 


P * (~Zl +Z2- ff(zi)) 
Zi — Z2 + Z 3 

—q * z 2 


(27) 


where g(z\) = mi * Z\ + 1/2 * (mo — mi) * (|zi + 1| — \zi — 1|), with the parameters taken values as p = 9.78, 
q = 14.97, too = —1.31 and mi = —0.75, which implies that the intrinsic node dynamics (without coupling terms) 
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have a double-scrolling chaotic attractor 35]. Let P = T = G = J 3 , where I 3 stands for the identity matrix of three 
dimensions. To estimate the parameter /3 in the QUAD condition, noting the Jacobin matrices of / is one of the 
following 


' -2.445 

9.78 

0 ' 


' 3.0318 

9.78 

0 ' 

1 

-1 

1 

,a 2 = 

1 

-1 

1 

0 

-14.97 

0 


0 

-14.97 

0 


then we estimate f}' = a — A max ((A 2 ) s ) = a — 9.1207, where 9.1207 is the largest eigenvalue of the symmetry parts 
of all Jacobin matrices of /. 

The possible coupling graph topologies are shown in Figj5] To select the pinned nodes, we add an extra virtual 
node (on behalf of s(t )) to the original network, which has a few links to the node that are pinned, then have an 
extended network. From the results in 0 , one can see that if every extended network topology among the switching 
is strongly connected and the duration time at each network topology is sufficiently long, then there always exist 
positive matrices P(-), coupling gain c and pinning gain e such that condition (1 1 21 ) holds, which implies that the 
system with persistent coupling and control can be stabilized. Hence, according to this viewpoint, one node from 
every strongly connected component is picked to be pinned. Noting topologies given in Fig [5] are all connected, every 
possible pinned node set is applicable. 

In this simulation, the graph topologies and pinned nodes of coupled system m switch among these four states, 
as shown in Fig[5](a)-(d) respectively, induced by a homogeneous Markov chain, o t . L(o t ) is picked as the Laplacian 
of the graph C?(cr t ), where each link has uniform weight 1. Here, we pick the state space of the Markov chain o t by 
{1,2, 3,4}, and its transition matrix is given by 


T = 


-10 

7 

0 

4 


6.5 

-10 

1 

6 


0 

3 

-10 

0 


3.5 

0 

9 

-10 


It can be seen from the transition matrix of a t , the expected sojourn time in each graph follows an exponential 
distribution with parameter 0.1. In the following, we pick a = 10 and /? = 0.8803. 

The ODEs (1112b are numerically solved by the Euler method with time step 0.001 (sec) and the time duration of 
the numerical simulations is [ 0 , 10 ] (sec). 


5.1. Continuous monitoring 

We give two examples to illustrate the updating rules (Il3l) and (IT4b respectively. Pick <5 = 0.03, c = 20, a = b = 
0.5 and e = 0.5. Then, it can be verified that {al m — kL{u) — keD(u)} s are negative definite and so the matrix 

inequality (fl2l > is satisfied. And we have = 0.2736. 

We employ rule (IT3l) . Fig. a shows the dynamics of each components of the 10 nodes and Fig. [ 8 ] shows the 
dynamics of V(t). All show that the coupled system i|T] [2} is stable. Similarly, we also employ the triggering event 
rule dTTb . Fig. [5]illustrates the dynamics of each components of all nodes, and Fig. [ 8 ]illustrates the dynamics of V ( t ). 
One can see that the coupled systems (1112b is asymptotically stable at certain chaotic homogeneous trajectory. 


5.2. Discrete monitoring 

In this subsection, we illustrate the discrete-time monitoring strategies as described in Theorem[2] In these exam¬ 
ples, we also take c = 20, a = b = 0.5, and e = 0.5. Fig| 6 ]shows the dynamics of each components of the 10 nodes 
and Fig| 8 ]shows the dynamics of V(t) under rule ( 125b in Theorem[2] All of them show that the coupled system (1112b 
is stable. By employing rule (l26l > in Theorem[2] Fig. [7]illustrates the dynamics of each components of all nodes and 
Fig. [ 8 ] illustrates the dynamics of V ( t ). These plots show that the coupled system (1112b is asymptotically stable. 

Furthermore, as shown in Figs |9(a)|9(b)| the events of updating the diffusion and pinning terms in the discrete 
monitoring strategy is much more than the continuous monitoring strategy, as we expected as Table Q] As a trade¬ 
off, the performance of the discrete monitoring in terms of convergence rate of V ( t ) is higher than the continuous 
monitoring, as shown by Fig | 8 ] Further, it can be seen from Fig. [ 8 ]that the rules ( 1 1 3b and (1251 ) have higher convergence 
rates than <U 1 and d26l > and close to the original coupled system with simultaneous diffusion and pinning, as a reward 
of high event frequency of updating diffusion and pinning terms, as shown in Figs|9(a)||9(b)| 
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(a) Pinned set {2, 6} 


(b) Pinned set {5, 8} 




(c) Pinned set {2, 6} 


(d) Pinned set {2, 5} 


Figure 3: The topologies of the graph of the coupled system and the pinned set. 


6. Strength, limitations of the work and orients of future research 


Event-triggered algorithm is a new issue in the coordination control. Despite attracting increasing interests re¬ 
cently, there are a small number of papers, for instance 0 - 0 , which were concerned with pinning control of 
networks with event-triggered algorithm. Moreover, all of them can not handle the scenario considered in this paper. 

12 M studied distributed event-triggered mechanism for pinning control of networks with static topology, while in our 
paper, we studied pinning networks with Markovian switching topologies and Markovian switching pinned node set. 
Il29ll30h investigated event-triggered pinning control of networks with static and switching topologies. However, from 
the sufficient conditions for complete synchronization that were given in E2l and 0, the pinned coupled system 
with each possible topology among the switching topologies should be able to stabilize the coupled system; in com¬ 
parison, in our work, there may exist some network topology and pinned node set in the state space of the Markovian 
chain that cannot stabilize the coupled subsystem. Moreover, they have not considered the switching of pinned node 
set, which were taken into consideration in the present paper. 

However, there are a few limitations of the present methods. First, the present study assumes there is no delay 
in information transitions. But real networks may have limit bandwidth limitation that will cause delays in message 
delivery. An interesting future research may take the time-delays into consideration. Second, in this work, we assume 
that the possible graph topologies and corresponding coupling matrices are already given and induced by a Markovian 
chain. It is sufficient for constructing a Lyapunov function to prove the stability of system. But in reality, the coupling 
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Figure 4: For continuous monitoring with event triggering ED, the dynamics of components of the coupled system ( 11121 >. 


weights of every possible graph topology may be time-varying, it is an important issue, and will be addressed in the 
near future. 

Third, this work mainly shows that if the linearly coupled system with persistent diffusion and control can be 
stabilized, then the proposed event-triggered rules can stabilize the system, too. It is important to extend the model 
to nonlinear cases. In 1361-1381. the authors proved that fuzzy systems are universal approximators for nonlinear 
dynamic systems. Hence, applying event-triggered strategies to fuzzy systems can be seen as a modest step. For this 
issue, we refer readers to [39]-[i4lj] . Recently, [42] proposed a centralized event-triggered communication scheme for 
networked Takagi-Sugeno fuzzy systems, while distributed event-triggered algorithms for fuzzy systems are absent. 
This also leads to interesting orients of our future work. 


7. Conclusions 

In this paper, event-triggered configurations and pinning control are employed to realize stability in linearly cou¬ 
pled dynamical systems with Markovian switching in both coupling matrix and pinned node set. Two monitoring 
scenarios are considered. For continuous monitoring, each node observes its neighborhood’s state and the target’s 
state (if it is pinned) in an instantaneous way to determine the next triggering event time for updating state informa¬ 
tion. Instead, for discrete monitoring, each node can only obtain the state information at the event time or switching 
time of the underlying Markov chain to predict the next triggering event time for updating state information. Once 
an event for a node is triggered, the diffusion coupling term and feedback control term of this node is updated. Event 
triggering criteria are derived for each node that can be computed in a parallel way. For both scenarios, it is proved 
that the coupled system can realize stability and the rule of piece-wise constant diffusion and pinning (if pinned) terms 
can efficiently reduce the computation load of the networked system, in comparison to the original coupled system. 
In addition, the discrete monitoring strategy can also reduce the communication load as well. Zeno behaviors can be 
proved excluded by proving the positivity of the lengths of the inter-event time intervals for some rules. Simulations 
are given to verify these theoretical results. 
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Figure 5: For continuous monitoring with event triggering 03- the dynamics of components of the coupled system 01121 ). 
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